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We study the nonequilibrium transport through a quantum dot coupled to normal and superconducting 
leads. We use the modified second-order perturbation theory to calculate the differential conductance and 
the local density of states at the quantum dot. In the strong but finite Coulomb interaction regime, the 
differential conductance shows an anomalous peak not at a zero bias voltage but at a finite bias voltage. 
We also observe an additional Kondo resonance besides the normal one in the local density of states, where 
the former is caused by nonequilibrium Andreev tunneling via the normal Kondo resonance. We explain 
that this specific Andreev tunneling gives rise to the anomalous peak in the differential conductance. Since 
the Andreev tunneling via the Kondo resonance is suppressed with increasing temperature, the anomalous 
peak in the differential conductance disappears at high temperatures. 
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In recent decades, the Kondo effect has received renewed 
interest in the context of mesoscopic systems fabricated with 
advanced nanotechnology. Particularly, the nonequilibrium 
transport through a quantum dot (QD)'' coupled to two nor- 
mal leads has been the subject of intense study, because 
interesting many-body effects are presumed to exist under 
nonequilibrium steady-state conditions. Furthermore, if such 
a QD is coupled to superconducting and normal leads (the 
N-QD-S system), more fascinating phenomena, such as the 
Andreev reflection, can be observed. 

In order to elucidate the interplay between the Andreev 
reflection and the Kondo effect, extensive experimental--^* 
and theoretical"*"'^' studies of N-QD-S systems have been 
conducted. Although some theoretical studies have already 
addressed the nonequilibrium transport in N-QD-S systems 
focusing on the influence of the Kondo effect,"'"'^' there re- 
main several important issues to be resolved. For example, 
there is still controversy about the zero-bias anomaly in the 
differential conductance: some studies have shown the ab- 
sence''^' (existence^' of a zero-bias anomaly in the Kondo 
limit. Moreover, most theoretical studies of the nonequilib- 
rium transport have been concerned with the limit of infinitely 
strong Coulomb interaction at the QD.**"^* However, since suf- 
ficiently large Coulomb interactions exclude the supercon- 
ducting correlations at the QD, no interplay between the An- 
dreev reflection and the Kondo effect may be observed under 
such an assumption. Thus, it is necessary to consider the in- 
termediate interaction regime. 

In order to clarify the interplay between the Andreev re- 
flection and the Kondo effect, we investigate the nonequi- 
librium transport in the N-QD-S system over a broad range 
of Coulomb interactions. For this purpose, we extend the 
modified second-order perturbation theory,"''^' which gives 
the self-energy interpolating between the weak and strong 
coupling limits, to a nonequilibrium case by employing the 
Keldysh technique. This method enables us to study the 
nonequilibrium system in the intermediate interaction regime. 
From the results of the differential conductance and the local 
density of states (LDOS) at the QD, we elucidate how the 
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Kondo effect affects the Andreev transport at a finite bias. 
In particular, a unique phenomenon in transport is observed 
in the intermediate regime: the nonequilibrium Andreev tun- 
neling via the normal Kondo resonance causes an anoma- 
lous peak in the differential conductance not at a zero bias 
voltage but at a finite bias voltage. Its origin is discussed in 
detail from the viewpoint of the competition/cooperation be- 
tween the Kondo and superconducting correlations. Note that 
since the anomalous peak appears only in the intermediate 
regime and seems to disappear in the limit of the infinitely 
strong Coulomb interaction, the results obtained here are con- 
sistent with most of the previous ones reported in this extreme 
limit.'*-'" 

We describe the N-QD-S system in terms of an Anderson 
impurity coupled to both the normal metal and the BCS su- 
perconductor, H - Hqd + //n + ^^s + Hj^ + Hjs, where 
^^QD = ^dTjo-dl-dg. + Udl^d^d^^d^ represents the Hamilto- 
nian for the QD. Hn = Zai^k - ^^N)c^fJ|^^c^|^^ and Hs = 

mal lead and superconducting lead, respectively (A is the su- 
perconducting gap). We assume that the bias voltage V is ap- 
plied only to the normal lead: - eV and /is = 0. The 
tunneling Hamiltonians between the QD and the two leads are 
^^TN = ZkAtNcli^^d^+H.c.) and//TS = ZkAtsclf,^d^-\-H.c.), 
where tN(S) is the tunneling amplitude between the QD and the 
normal (superconducting) lead. We consider the wide-band 
limit of electrons in the leads, in which Fncs) = ^fli(S)Pi<(S,) 
fined at A = represents the resonance strength between the 
QD and the normal (superconducting) lead. 

The time-dependent current flowing from the normal (su- 
perconducting) lead into the QD, /jv(S), is described by 

lN(S}(t) = (ie/h) Yjkcr tN{S)cl;^^s)ka-^'^^'r(') + ^-C Sincc we con- 
sider the nonequilibrium steady-state transport, the expecta- 
tion values of the current {In{S)} is time-independent. Using 
the standard Keldysh Green's function formalism, {In(S)) can 
be expressed in terms of the local retarded and lesser Green's 
functions at the QD, G^(cj) and G''{oj)^ i"*) The bold font rep- 
resents a 2x2 Nambu matrix. These Green's functions include 
the self-energies due to the Coulomb interaction U in addition 
to those due to the coupling between the QD and the leads. In 
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order to obtain the self-energies for U, we use the modified 
second-order perturbation theory"' in the Keldysh Green's 
function formalism which is a generalization of Cuevas's ap- 
proach"* in the equilibrium N-QD-S system to the system 
under nonequilibrium conditions. Following ref. 11, the mod- 
ified retarded self-energy j^^^ is given by a functional of the 

bare second-order self-energy such that S^^,,^ has 

an appropriate form in both the atomic limit, Fai s / f/ — » 0, 
and weak-interaction limit, U ITm,s 0. In a similar manner, 
we obtain the modified lesser self-energy, S^,^^."' 

We can calculate the bare second-order self-energies 
and £y2«d using the perturbation theory with respect to U in 
the standard Keldysh Green's function formalism. The only 
point to note in this approach is that the Hamiltonian for the 
QD is replaced by the following modified one when we cal- 
culate the self-energy: - Id Ho- dlda- + (A^^fj^f j -i- H.c.) + 
A{In + Is ), where the effective parameters ed and Ad are intro- 
duced to ensure that the self-energies have the correct forms 
in the atomic limit."' '■'^ In addition, we introduce the source 
term A coupled to the current operator in order to conserve the 
current, because the simple application of the second-order 
self-energy in the N-QD-S system without the source term 
may break the current conservation law: {1^} + {Is) + 0.'^' 
The effective parameters Aj, and ~A are then fixed under 
the following conditions, which respectively indicate the con- 
sistency of the electron number, the effective counter gap to 
the proximity effect at the QD, and the current conservation: 
<«> = <«>o, A^ = [i:'f;(0)]i2 and {U) + </s> = 0. In the first 
equality, («) and («)o are the expectation values of the elec- 
tron number at the QD with the full Hamiltonian and with the 
modified one. 

First, we discuss how the differential conductance dZ/dV 
is affected by the Coulomb interaction. Figure 1(a) shows the 
differential conductance as a function of the bias voltage Y for 
several values of IJ with the symmetric coupling F^ /Fai - 1 . 
We only show dZ/dV for a positive bias voltage since it is an 
even function of W in this case (ej + V 12 - [is - 0). For 

= 0, the zero-bias conductance takes its maximum value 
Ae~ jh, where the Andreev tunneling between the QD and the 
superconducting lead is balanced by the normal electron tun- 
neling between the QD and the normal lead. This balance is 
ensured by the condition of Ts ITn - 1 with U - Q. With 
the enhancement of the Coulomb interaction, the differential 
conductance decreases gradually in the entire voltage range, 
because the Coulomb interaction suppresses both the Andreev 
tunneling and the normal electron tunneling, but in different 
ways: the interaction makes the former much weaker than the 
latter This is evident in renormalized parameters, f and Fj , 
which are defined for V - Q and T = as 



fs =z(Fs +[Re2:j/(w)]i2), 



(1) 



where z is the renormalization factor z — (1 + F^/A - 
d |^ReIl^(w)jjj /dw]^ q)"'. Figure 1(b) shows the U depen- 
dence of these renormalized parameters. It is seen that 
decreases more rapidly than f with an increase in U. For 
fs < Fiv, the coupling between the QD and the supercon- 
ducting lead effectively becomes weaker than the coupling be- 
tween the QD and the normal lead, making the Kondo singlet 
dominant at the QD. As a result, the system loses the above 
subtle balance in tunneling processes, leading to a global re- 
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Fig. 1. (Color online) (a) Differential conductance as a function of bias 
voltage V for several values of the Coulomb interaction U: £d/U = —0.5, 
A/Yf/ = 1, Tg /Fa; = 1, and ksT/Yfj = 0.01. (b) Renormalized parameters 
Fa/ and Tg as a function of the Coulomb interaction U: £j/U = —0.5, 
A/Fw = 1, Fs/Fw = 1, sV/Tn = 0, andArgr/Fjv = 0. 



duction in conductance. Note here that the differential con- 
ductance markedly changes its voltage dependence with an 
increase in U (Fig. 1(a)). In particular, the monotonic curve 
of the differential conductance for U/Tn - gradually de- 
velops a double -peak structure, as seen for U/Fn - 10. As 
U further increases, the two peaks become more prominent. 
The peak in the vicinity of V = becomes sharper and 
gradually approaches a zero-bias voltage: e.g., it is located at 
eV/Ff^ ^ 0.05 for f//F/v - 20. The anomalous sharp peak 
is commonly observed when the Kondo effect is dominant 
at the QD (Ts/^n < 1 and {//Fa, » 1) though its height 
decreases with increasing U. On the other hand, the broad 
peak at eV/TN - 0.6 is insensitive to the increase in U. The 
above novel characteristic feature is one of the main findings 
of this study, which is quite different from an ordinary zero- 
bias anomaly in a QD coupled to two normal leads. Note that 
with increasing U, the anomalous peak gets close to the zero- 
bias voltage but decreases in height. Eventually, the anoma- 
lous peak seems to disappear in the limit of f/ —> oo, which is 
consistent with the absence of the zero-bias anomaly in refs. 
4, 5, 6, and 10 where an infinitely strong Coulomb interaction 
is assumed. The monotonic decrease in the zero-bias conduc- 
tance with increasing U is also consistent with the results of 
equilibrium studies. "• 

We also show the differential conductance in the cases with 
asymmetric couplings (Fj /r^ + 1) and a strong Coulomb in- 
teraction [//Fa, = 20 in Fig. 2(a). For F^/Fa- > 1 (< 1), the 
sharp peak near the zero-bias voltage increases (decreases) 
in height and changes its position toward a slightly higher 
(lower) bias voltage. To be more specific, let us denote the po- 
sition and height of the sharp peak as Va and Ga, respectively, 
which are plotted in Fig. 2(b) together with the renormal- 
ized parameters. A significant feature is that the peak position 
(cVa) is approximately given by Fa, for small F^'s, where Fa, 
corresponds roughly to the Kondo temperature. Also note that 
the peak height Ga approaches zero for Fj ^ because the 
Andreev tunneling completely disappears in that limit. These 
facts indicate that the origin of this characteristic peak is re- 
lated to both the Kondo effect and the proximity effect. Note 
that the other broad peak is also affected by the change in Fj . 
For smaller Fj /Fa, values, the broad peak at eV/FA, - 0.6 in 
the symmetric case (Fj /Fa, = 1) shifts toward the gap edge. 
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Fig. 2. (Color online) (a) Differential conductance as a function of bias 
voltage y for several values of Y$ and U /Tn = 20. The other parameters 
are the same as those in Fig. 1(a). (b) Anomalous sharp peak position Va 
and its height Ga as a function of Ts . The parameters are the same as those 
in (a). Here, Go = I6e^ /h. The renoiTnalized pai'ameters defined for V = 
and r = 0, Fiv and Ts , are also shown. 

eV/r^ = 1. In contrast, for larger Fs /Fat values, the broad 
peak shifts in the opposite direction and merges with the sharp 
peak near the zero-bias voltage. All the above characteristic 
features of the conductance will be clearly explained below 
in terms of the interplay between the Kondo/proximity effects 
under nonequilibrium conditions. 

Let us now discuss how the LDOS at the QD, p(w), in the 
strong Coulomb interaction regime changes its profile as bias 
voltage increases. The computed results are shown in Fig. 3. 
Let us first look at the LDOS in the equilibrium state {V = 0) 
for U/Fn - 20, which is compared with that in the case of 
U/Tn = shown in the inset. The two resonances emerging 
inside the superconducting gap for U/r^ - represent the 
Andreev resonant states, which are approximately described 
by the superposition of two Lorentzians with renormalized 
parameters: they are located at w = with a width Ff^. It 
is seen that the Andreev resonance for U/T;^ = is strongly 
renormalized and merged into a single sharp Kondo resonance 
at the Fermi level for U/Fn = 20 because f s decreases more 
rapidly than (see Fig. 1(b)). The resulting Kondo reso- 
nance is sensitive to the change in bias voltage, so that the 
LDOS for U/Tn = 20 significantly alters its shape with an 
increase in V. Indeed, as bias voltage increases, the position 
of the Kondo resonance shifts following the Fermi level of 
the normal lead fiN- e.g., for iun/^n = eV/r^ = 0.2, the peak 
position is located at w/F/y ^ 0.2. This implies that the ordi- 
nary Kondo resonance is formed by conduction electrons in 
the normal lead. A noteworthy feature for eV/F^i = 0.2 is that 
another resonance develops, though small, near the counter- 
position of the Kondo resonance a> ^ -fJ-N- This additional 
Kondo resonance was previously observed by Sun et al.,'^^ 
but has not been discussed in detail, particularly as regards 
its physical relevance to the transport properties. We will ad- 
dress this issue below and demonstrate that it indeed provides 
a source of the marked change in nonequilibrium transport 
properties. For eV/r^ = 0.6, the peak value of the addi- 
tional Kondo resonance is slightly increased while the nor- 
mal Kondo resonance is suppressed, leading to a broad two- 
peak structure, which is analogous to that in the U /F^v = 
case. The two-peak structure therefore indicates that the su- 
perconducting pairing state, which is strongly suppressed by 



Fig. 3. (Color online) LDOS at the QD for several values of bias voltage V. 
We set ej/U = -0.5, C//Fjv = 20, Fs/F^ = 1, A/Fa- = 1, and ksT/rN = 
0.01. The LDOS for eV/Tfj = is shown in the inset in comparison with 
the one for the noninteracting case (thin solid line). 

the large Coulomb interaction at V = 0, is revived and be- 
comes dominant over the Kondo singlet state. With further 
increase in bias voltage (eV/FN = 1,2), however, the An- 
dreev resonance is suppressed again since the applied bias 
voltage also makes the superconducting pairing state unsta- 
ble. In these cases, the weight of the LDOS is transferred to 
the region at w ^ +U/2, corresponding to the characteris- 
tic energy of charge excitations at the QD (data not shown), 
which indicates that the QD is in the local-moment (free-spin) 
regime. Therefore, we conclude that as bias voltage increases, 
the LDOS first exhibits a crossover from the Kondo-dominant 
regime to the superconducting-dominant regime, and then to 
the local-moment regime. 

The above discussions on the LDOS enable us to clarify 
the origin of the two peaks in the differential conductance in 
the large-?/ regime in Figs. 1(a) and 2(a). We first recall that 
the zero-bias conductance is substantially suppressed by U, as 
has been shown in Fig. 1(a). However, for a small bias volt- 
age, the effective cotunneling process is enhanced owing to 
the nonequilibrium Andreev tunneling via the Kondo reso- 
nance, which gives rise to the anomalous increase in the dif- 
ferential conductance at a small but finite bias voltage. Conse- 
quently, we can say that the anomalous increase found in the 
conductance near zero bias is a unique nonequilibrium phe- 
nomenon induced by the interplay of the Kondo effect and 
superconducting proximity effect. Here, we determine why 
the anomalous peak is located at eV/FN - 1. For f/ = in 
Fig. 1(a), the zero-bias conductance has a maximum value for 
f s /Fn = 1 • Namely, the maximum transport probability is re- 
alized when the width and position of the Andreev resonance 
have the same value. This specific correspondence holds for 
the zero-bias conductance even in the interacting case be- 
cause the low-energy physics of the N-QD-S system can be 
described by the local Fermi liquid theory." '-' For instance, 
as Fs becomes smaller than F^ with increasing U, the system 
loses the above specific balance, resulting in the suppression 
of the zero-bias conductance. Even under such conditions as 
Fs /Fn 1, however, the peak position of the Andreev reso- 
nance is shifted and determined by the bias voltage V instead 
of by Fs , as shown in Fig. 3. Accordingly, we can state the fol- 
lowing by generalizing the above correspondence to the non- 
linear case: In nonequilibrium cases with eV/F^ - 1, where 
the position of the Andreev resonance is the same as its width. 
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Fig. 4. (Color online) (a) Differential conductance and (b) LDOS as a func- 
tion of bias voltage V for several values of temperature T: UIT^ = 20, 
ejlU = -0.5, Fs/Fn = 1, and A/Fjv = 1. 



the transmission probability is enhanced, which leads to the 
anomalous peak of the differential conductance. In contrast, 
at eV/FN - 0.6, the superconducting pairing state is domi- 
nant and thus the Andreev reflection is enhanced, leading to a 
broad peak in the bias voltage dependence of the conductance. 
Namely, the two-peak structure found for the conductance in 
Fig. 1(a) clearly characterizes a crossover from the Kondo- 
dominant regime to the superconducting-dominant regime, 
and then to the local-moment regime. Even in the asymmetric- 
coupling cases shown in Fig. 2(a), we can see the increase in 
conductance at a small finite bias voltage, the origin of which 
is the same as that in the symmetric case. On the other hand, 
the broad peak at higher voltages depends on the asymmetry 
of the couplings. For F^ /F/^ = 0.5 (=2.0), the Kondo correla- 
tion is stronger (weaker) than that in the symmetric coupling 
case, so that a larger (smaller) voltage is necessary to make 
the superconducting correlations dominant, which shifts the 
peak position to a higher (lower) voltage side in Fig. 2(a). 
The above arguments clearly explain the origin of our new 
findings on conductance in this paper 

Finally, we show the temperature dependences of the dif- 
ferential conductance and LDOS. Figure 4(a) shows differ- 
ential conductance as a function of bias voltage for several 
temperatures T. It is seen that with increasing T, the sharp 
peak near eV/F^ = 0.05 begins to collapse at fc^r/F/v ^ 0.05, 
and completely disappears at ksT/Tj^ = 1. This is because 
the Andreev-normal cotunneling process due to the Kondo ef- 
fect is suppressed at high temperatures. In order to confirm it, 
we also calculate the temperature dependence of the LDOS 
for eV/FN = 0.05, as shown in Fig. 4(b). We see that the 
Kondo resonance begins to collapse at IVm ^ 0.05, and 
its height for ksT IT^ = 0.1 becomes smaller than half of 
that in the low-temperature limit. Note that at this character- 
istic temperature ksT lY^ - 0.05, the conductance peak at 
cVITn - 0.05 is also suppressed, ensuring that the anomalous 
peak in the conductance is indeed caused by the Kondo effect. 
These characteristic eV and kgT correspond approximately to 
Tn- In contrast, the broad hump at eVIYt^ - 0.6 is less sensi- 
tive to the temperature because it is caused by the proximity 
effect but not by the Kondo effect. Further increase in temper- 
ature finally suppresses both the Kondo effect and the proxim- 
ity effect, so that there is no resonance in the LDOS for both 



eVIYjsi — 0.05 and 0.6 at a high temperature of kgT jYf] - I. 
This naturally explains the decrease in the differential con- 
ductance in the entire V range at higher temperatures shown 
in Fig. 4(a). 

In summary, we have investigated the nonequilibrium 
transport in a QD system coupled to normal and supercon- 
ducting leads where the Andreev reflection plays an impor- 
tant role. Using the modified second-order perturbation the- 
ory, we have elucidated that the differential conductance de- 
velops two characteristic peaks in its bias voltage dependence 
as the Coulomb interaction is enhanced. In particular, the 
peak near the zero-bias voltage characterizes a unique phe- 
nomenon caused by a nonequilibrium Andreev tunneling via 
the Kondo resonance. It has indeed been confirmed that the 
sharp peak disappears with increasing temperature, reflect- 
ing the suppression of the Kondo resonance. The two-peak 
structure found for the voltage dependence of the conductance 
characterizes a crossover from the Kondo-dominant regime to 
the superconducting-dominant regime, and then to the local- 
moment regime. In real experiments, different patterns in the 
temperature dependence of the peaks may be observed, which 
can give strong evidence of the interplay between the Andreev 
reflection and the Kondo effect at a finite bias. We think that 
the transport experiments proposed here will provide an im- 
portant test bed for a deeper understanding of the nonequilib- 
rium Kondo/ Andreev physics in QD systems. 
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